† Corresponding author. E-mail:
Using the first-principles method based on the density functional theory (DFT), the structures and electronic properties of different gas hydrates (CO2, CO, CH4, and H2) are investigated within the generalized gradient approximation. The structural stability of methane hydrate is studied in this paper. The results show that the carbon dioxide hydrate is more stable than the other three gas hydrates and its binding energy is −2.36 eV, and that the hydrogen hydrate is less stable and the binding energy is −0.36 eV. Water cages experience repulsion from inner gas molecules, which makes the hydrate structure more stable. Comparing the electronic properties of two kinds of water cages, the energy region of the hydrate with methane is low and the peak is close to the left, indicating that the existence of methane increases the stability of the hydrate structure. Comparing the methane molecule in water cages and a single methane molecule, the energy of electron distribution area of the former is low, showing that the filling of methane enhances the stability of hydrate structure.
Gas hydrate is a kind of clathrate that is formed by water and gas molecules under certain pressure, temperature, salinity and chemical composition conditions.[1] The natural gas hydrate reservoirs are outstanding for their huge reserve, extensive distribution, high energy density and high caloric value.[2] The energy efficiency of the gas hydrate is ten times that of coal.[3] It is thought that it is going to be an important potential source of energy in the future.
With the development and innovation of molecular simulation technology, the structures and properties of gas hydrate have been extensively investigated by using molecular dynamics (MD), Monte Carlo (MC) and density functional theory (DFT) methods.[4–6] Duan et al.[7] chose the constant-pressure and constant-temperature (NPT) ensemble, the consistent valence force field and the Ewald method for molecular dynamics simulation of sI methane hydrate supercell system. They also analyzed the changes of system conformation, radial distribution function, system potential energy and mean square displacement under different pressure conditions have been analyzed. Tang et al.[8] employed the DFT-D method to study the phenomenon of water cage merging in the process of methane hydrate nucleation. They analyzed the structures, stabilities and Raman spectra of double cages and multi cages under the lowest energy. Song et al.[9] calculated the initial structure of methane hydrate without a single or two water molecules and optimized the unrestricted structure at a B3LYP/6-31 + G(d, p) level. The BSSE method was used to calculate the interaction energy between water and methane. They also studied the relationship between the cage structure and stability of methane hydrate without water molecules. Cao et al.[10] systematically calculated the stabilities of 18 alkane guest molecules in two kinds of water cages (51262 and 51264) by density functional theory at a B97-D/6-311 ++ G(2d, 2p) level. They simulated Raman spectra of nine kinds of alkane hydrates to study the stability between the large guest alkane molecule and the water cage. They also studied the change of the vibrational frequency of the guest molecule in the hydrate with the object molecular category. According to the first principles based on the density functional theory method, Xu[11] analyzed the interaction between the guest molecules and water cage architecture. The binding energies of methane hydrate and ethane hydrate under the structure of sI and sII were calculated respectively. They also studied the stability of each system and the influence of tetrahydrofuran (THF) on the system structure.
Through the above investigation, we could see that the structural properties of methane hydrate and the change of the system energy have been calculated. Our work is based on these previous researches. In this paper, we calculate the binding energies of several common gas hydrates, and also analyze the methane hydrate structural stability from the view of density of states(DOS). In this paper, we not only calculate the partial electronic density of states (PDOS) of a water cage with/without the inner methane molecule, but also analyze the PDOS of single methane and methane in a water cage, to further verify the structural stability of methane hydrate. Among all kinds of gas hydrates, methane hydrate is outstanding for the huge reserve and extensive distribution. Therefore, we choose methane hydrate as the main research object in this paper. Through this research we hope that the mechanism can be revealed for quantum mechanics of gas hydrate structure stability, especially methane hydrate, and the theoretical guidance can be provided for the research and utilization of gas hydrate.
Based on the first principle method of DFT,[12,13] the geometric structure and electronic properties of different gas hydrates are calculated by the quantum mechanical program CASTEP which is developed by the condensed matter physics research group of the University of Cambridge. The program reliability has been verified by a large number of practical calculations. The calculation of the CASTEP program is based on the total energy pseudo potential method and the plane wave expansion,[14] which requires the system to be periodic.
In this paper, the initial structure of gas hydrate is derived from the x-ray diffraction experiment.[15] The space group of methane hydrate is PM3N, and the lattice constant is a = b = c = 1.162 nm. The calculation is carried out in the reciprocal space, so the volume of real space needs to be converted into the inverted lattice space. The exchange correlation function is based on the generalized gradient approximation (GGA-PBE).[16] The K point[17] grid size is 2×2×2. Self-consistent field cycle converges to 1.0×10−5 eV. The BFGS[18] optimization algorithm is used to optimize the geometry. After geometry optimization of the gas hydrate structures, the parameters of the crystal structure at ground state are obtained. The total energy (ET) of methane hydrate is −21275.69 eV. The forces acting on each atom are less than 0.1 eV/nm. The value of plane wave cut-off energy is 100 eV. The convergence test is carried out by changing the cut-off energy. The results show that these settings are sufficient to ensure the accuracy of the calculation. The structures of gas hydrate with/without an inner methane molecule are shown in Figs.
In this paper, we use 1×1×1 periodic super cells as the initial structure. The study is divided into two cases, one is the water cage with an inner methane molecule (Fig.
In order to quantitatively reflect the difference between the two structures, the distances from the cage center to the oxygen atoms of the water cage are calculated (Table
To analyze the stability of gas hydrate, we calculate the binding energies of different gas hydrates. The formula is as follows:
Figure
Then we study the difference between water cages with/without the inner methane molecule from the view of DOS. As shown in Fig.
On the one hand, we compare the PDOSs of the water cage with/without the inner methane. As shown in Fig.
On the other hand, we conduct a comparison of PDOS between single methane and methane in a water cage (Fig.
Using the first principles method based on density functional theory, we have studied the structural stabilities of different gas hydrates, the binding energy, ground-state structures and electronic properties of different gas hydrates. A comparison of the two kinds of structures shows that water cages experience the repulsion from inner gas molecules, which makes the hydrate structure more stable. The analyses of the electronic properties of the two water cages show that the energy region of hydrate with inner methane is lower, and the peaks are closer to the left, indicating that the existence of methane increases the stability of the hydrate structure. The analysis of the density of states shows that when the methane exists, the distribution of s/p electrons of the water cage shifts to the left and the energy is reduced. In addition, the electron distribution of methane shifts to the left and the energy decreases after binding to the water cage, which further illustrates that hydrate structure is stable in the presence of methane. The binding energy of methane hydrate is −0.58 eV, which is higher than that of carbon dioxide hydrate (−2.36 eV) and lower than that of hydrogen hydrate (−0.36 eV), showing that the structure stability of carbon dioxide hydrate is higher than that of methane hydrate, while hydrogen hydrate is the most unstable among them.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] | |
[15] | |
[16] | |
[17] | |
[18] |